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where  a  hazard  would  exist  as  well  as  the  maximum  distance  down  wind  that  a 
hazard  would  exist.  Depending  upon  the  input  data  supplied,  the  program  will 
also  predict  the  amount  of  water  that  would  be  required  for  sufficient  dilu- 
lion  to  eliminate  the  hazard,  the  effects  of  changing  variables  (as  for 
example,  the  effects  of  changes  in  wind  speed,  or  in  air  temperature  or  in 
atmospheric  stability),  and  the  effect  of  the  passage  of  time  (up  to  twenty- 
four  hours)  if  weather  conditions  are  updated  hourly.  In  the  latter  case, 
allowances  are  maJe  for  composition  changes  in  the  remaining  pool  because  of 
preferential  evaporation  of  components,  possible  absorption  of  water  from  the 
atmosphere,  and  absorption  of  (and  reaction  with)  carbon  dioxide  from  the 
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EVAPORATION  AND  DISPERSION  OF  HAZARDOUS  MATERIALS 


Introduction: 

This  work  was  undertaken  In  order  to  develop  a  state-of-the-art  predictive 
•ode  1  for  the  evaporation  and  dispersion  of  hazardous  liquids,  especially  hydro- 
sine  family  fuels  and  propellants,  which  sight  be  accidentally  spilled. 

The  Importance  of  the  problem  has  been  adequately  shown  previously  (1,  2). 

Some  earlier  work  by  Ills  and  Springer  (3)  approached  the  problem  from  the  stand¬ 
point  of  pure  materials.  However,  the  materials  of  particular  Interest  here  are 
not  always  transported  or  used  as  pure  materials,  hence  In  this  work  it  has  been 
necessary  to  Incorporate  modifications  of  the  earlier  work  (3)  to  Include  mixture 
behavior . 

An  additional  consideration  has  been  the  incorporation  of  model  elements  re¬ 
lating  to  the  reactions  of  the  hazardous  materials  with  the  environment,  especially 
water  absorption  from  the  air,  carbon  dioxide  absorption  from  the  air  and  decay  of 
the  evaporated  material  during  Its  transport  down  wind. 

Since  the  users  of  this  predictive  model  will  ordinarily  wish  to  define  the 
down-wind  area  wherein  a  hazard  would  exist  in  the  case  of  a  spill,  the  predictive 
scheme  includes  down  wind  dispersion  and  hazard  corridor  definition. 

All  of  the  several  model  elements  have  been  assembled  and  combined  In  a  FORTRAN 
computer  program  which  will  compute  evaporation  rates  and  limits  of  the  hazardous 
region,  that  Is,  the  maximum  distance  away  from  the  centerline  where  a  hazard  would 
•mist  as  well  as  the  maximum  distance  down  wind  that  a  hazard  would  exist.  Depend¬ 
ing  upon  the  input  data  supplied,  the  program  will  also  predict  the  amount  of  water 
that  would  be  required  for  sufficient  dilution  to  eliminate  the  hazard,  the  effects 
of  changing  variables  (as  for  example,  the  effects  of  changes  in  wind  speed,  or  in 
air  temperature  or  in  atmospheric  stability),  and  the  effect  of  the  passage  of  time 


1 


2 


(up  Co  twenty-four  hours)  if  weather  conditions  are  updated  hourly.  In  the  latter 
case,  allowances  are  made  for  composition  changes  In  the  remaining  pool  because  ot 
preferential  evaporation  of  components,  possible  absorption  of  water  from  the  atmos¬ 
phere,  and  absorption  of  (and  reaction  with)  carbon  dioxide  from  the  atmosphere. 

A  number  of  computations  were  made  In  order  to  determine  the  effects  of  changes 
In  the  several  input  parameters.  The  results  of  these  are  presented  in  the  form 
of  correction  factors  which  can  be  used  to  correct  from  certain  "base  cases",  and 
thus  provide  for  rapid  field  estimation  when  needed. 
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Nature  of  the  Hazardous  Material*  Considered: 

There  ere  three  hazardous  (toxic)  materials  of  interest  in  this  work;  they 
being:  hydrazine,  monomethyl  hydrazine  (MMH) ,  HjNNHCHj;  and  unsymetr leal 

dimethyl  hydrazine  (UDMH) ,  ^NNfCHj),.  Of  these  materials,  as  used  by  or  transported 
by  or  for  the  U.S.  Air  Force  only  the  MMH  is  In  pure  form.  The  UDMH  Is  used  as  an 
equal  parts  (by  mass)  mixture  with  hydrazine.  This  mixture  Is  sometimes  referred 
to  as  "Aerozlne  50".  Hydrazine  is  also  used  in  a  water  solution  of  70  percent  (by 
mass)  hydrazine. 

The  vapor  phase  toxicity  of  these  materials  has  been  defined  (4),  and  the  Short 
Term  Public  Exposure  limits  (STPL)  are  given  as  7.0  mg/m-*,  2.8  mg/m*  and  38.2  mg/m* 
for  hydrazine,  MMH  and  UDMH  respectively. 

Hydrazine  and  both  of  the  derivatives  are  miscible  in  all  proportions  with 
water.  Not  only  are  they  miscible,  they  arc  extremely  non- ideal  and  fit  the  defini¬ 
tion  of  "hygroscopic",  that  is,  they  very  readily  take  up  water,  and  the  equilibrium 
partial  pressure  of  water  vapor  over  the  solutions  is  much  lower  than  that  predicted 
by  Raoult's  law.  Some  limited  vapor-llquld  equilibrium  data  have  been  published 

(5,  6),  and  neither  the  hydrazine  and  water  system  nor  the  UDMH-water  system  obey 

% 

the  Marqules  equation  (7,  8).  This  non-ldeallty  is  with  respect  to  the  liquid  solu¬ 
tions.  There  are  no  indications  of  any  particular  non-ldcality  of  the  vapor  mixtures. 

The  three  hydrazine  species  each  have  somewhat  different  physical  properties, 
with,  in  most  cases,  the  properties  of  MMH  being  intermediate  between  those  of  hy¬ 
drazine  and  UDMH.  Most  of  the  relevant  physical  properties  of  the  pure  liquids 
have  been  collected  into  a  single  publication  (9). 
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U  “  viscosity 
p  “  density 

D  *  diffusivity  of  diffusing  component  in  air 
The  appropriate  properties  for  the  Schmidt  number  evaluation  is  the  gas  film  adjacent 
to  the  liquid  surface. 

The  evaporative  flux,  then  is  given  by: 

Qm  *  k«<Cl  -  CA)  (2) 

where: 

“  evaporative  flux 

Cj  •  equilibrium  concentration  of  diffusing  material  immediately  adjacent 
to  the  liquid  surface 


Ca  “  concentration  of  diffusing  component  in  the  bulk  air  stream 
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For  any  concentration,  C  -  p/KT  If  the  ideal  gas  low  holds,  where: 
p  ■  partial  pressure  of  the  component 
R  ■  universal  gas  constant 
T  *  temperature 

It  was  assumed  that  each  transfering  component  will  transfer  in  accordance  with 
equation  (2),  since  it  is,  essentially,  a  definition  of  the  mass  transfer  coefficient. 
For  an  evaparating  hydrazine,  CA,  that  is,  the  bulk  concentration  in  the  air,  is 
zero  at  all  times,  and  that  the  interface  composition  may  be  determined  by  assum¬ 
ing  it  to  be  an  equilibrium  concentration.  For  water,  which  may  transfer  in  either 
direction,  the  bulk  concentration  varies  with  the  humidity,  while  the  interlace  com¬ 
position  will  be  an  equilibrium  value.  For  carbon  dioxide,  it  was  assumed  that  the 
composition  was  0.03  percent  in  the  air  and  was  zero  at  the  interface. 

The  Schmidt  numbers  as  well  as  the  equilibrium  concentrations  are  functions  of 
temperature,  with  the  latter  being  very  strong  functions  for  each  component  (except 
carbon  dioxide).  The  estimation  methods  are  presented  in  Appendix  A. 

The  temperature  immediately  adjacent  to  the  liquid  surface  was  assumed  to  be 
equal  to  the  liquid  pool  temperature,  while  the  mean  temperature  of  the  film  (for 
Schmidt  number  evaluation)  was  assumed  to  be  the  mean  of  the  liquid  pool  temperature 
and  the  bulk  air  temperature 

The  liquid  pool  temperature  depends  upon  the  relative  rates  of  heat  transfer 
to  and  from  the  pool.  The  transfer  to  the  pool  being  radiant  heating  from  sunshine, 
qt,  radiation  from  the  atmosphere,  qa,  convective  heating  from  the  air,  qh  and  con¬ 
vective  heating  from  the  ground  surface,  qg.  It  is  possible  that  qh  and/or  qg  could 
be  negative.  The  heat  losses  are  rcradiatlon  from  the  pool  to  space  qr ,  and  the 
latent  heat  of  vaporizing  mat. ’rials,  q^. 

The  latent  heat  term  is  juite  large,  and  is  itself  a  strong  temperature  func¬ 
tion  since  the  evaporation  rate  is  a  strong  taaperature  function.  Therefore,  it 
is  reasonable  to  assume  that  the  actual  pool  temperature  is  that  which  is  in  equlli- 
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brium  with  the  surroundings,  that  is,  the  rate  of  heat  gain  is  equal  to  the 
rate  of  heat  loss. 

Thus : 

Is  ♦  9a  ♦  <lh  +  9g  "  9r  +  ^1  O) 

Except  for  q„  and  qa,  each  term  In  equation  (3)  is  a  function  of  the  pool 
temperature,  and  so  it  is  possible  to  deduce  the  pool  temperature  which  would 
be  the  equilibrium  temperature  by  finding  that  temperature  which  will  satisfy 
equation  (3). 

The  evaluation  of  q8  was  done  by  an  algorithm  devised  to  determine  the 
solar  Insolation  rate  on  horizontal  surfaces  as  a  function  of  latitude,  longi¬ 
tude,  date,  time,  and  cloud  cover  (11).  The  algorithm  is  rather  lengthy  and 
is  not  reproduced  here.  The  absorptivity  of  the  pool,  Up,  was  assumed  to  be  0.96 
which  compaies  well  with  water  and  a  number  of  organic  liquids  (10,  12). 

The  atmospheric  long  wave  length  radiation,  qa,  was  evaluated  as  a  function 
of  the  temperature  of  the  air  and  the  atmospheric  emissivity,  F.a ,  as  given  by 
Sax  (13). 

The  convective  transfer  from  (or  to)  the  air  was  modelled  after  MacKay  (10) 
which  Is  a  heat  transfer  analog  of  the  same  mass  transfer  coefficient  given  earlier 
as  equation  (1).  The  form  for  the  analogous  equation  is: 

ha  "  kmpCp(Sc/Pr)0'67 

where : 

p  "  density  of  the  "film" 

Cp  •  heat  capacity  of  the  "film" 

Pr  ■  Prandtlc  no.,  Cpp/k  of  the  "film" 
k  -  thermal  conductivity 

The  film  properties  were  evaluated  at  the  mean  temperature  and  the  mean  composi¬ 
tion  of  the  interface  (equilibrium  composition  at  the  pool  temperature)  and  the 


bulk  air  condition. 
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l'he  convective  transfer  between  the  ground  and  the  pool  was  treated  accord¬ 
ing  to  Schmidt  and  Sllveston  (14).  This  Is  the  same  algorithm  used  by  11 le  and 
Springer  (3),  although  the  liquid  mixture  properties  had  to  be  evaluated  as  com¬ 
position  functions  as  well  as  temperature  functions.  These  properties  estimations 
and  others  are  presented  In  Appendix  A. 

The  re-radiation  from  the  pool  is  that  of  a  radiator  radiating  to  space.  The 
liquid  pool  cmissivity,  Ep,  is  assumed  to  be  0.96  (12). 

Finally,  the  latent  heat  of  vaporization  is  the  specific  latent  heat  of  each 
component  times  its  evaporative  flux.  When  the  several  terms  of  Equation  (3)  are 
rewritten  as  temperature  functions,  it  becomes: 

qs^oE.  ♦  ha(Ta  -  Tp)  ♦  hK(Tg  -  y  -  TpO  ♦  LX^  (5) 

where: 

C  "  Steffan-Boltzmann  Constant 
Each  term  of  this  equation  is  an  energy  flux. 

Equation  (5)  may  be  solved  for  the  pool  temperature,  Tp,  if  all  the  other 
terms  arc  known.  The  solution  requires  an  iterative  method,  but  is  readily  found 
by  the  N'cwton-Raphson  method. 

The  mass  transfer  flux,  as  given  by  Equation  (2),  is  seen  to  be  a  function  of 
the  equilibrium  concentration  of  transfering  component.  In  Equation  (2)  this  concen¬ 
tration  is  identified  as  Cj.  For  any  component,  i,  this  equilibrium  composition  is 
given  by: 

(Cl>i  *  *A/RTp  (6) 

where: 

^ "  activity  coefficient 
Pv  "  vapor  pressure 

x  -  mole  fraction  of  component  in  the  liquid 
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Thu  activity  coefficient  is  predicted  by  methods  indicated  in  Appendix  A. 
Equation  (6)  watt  not  applied  to  carbon  dioxide,  but  rather,  the  equilibrium  compo¬ 
sition  was  a«t  at  zero. 

The  bulk  phase  composition,  CA,  was  set  at  zero  for  the  hydrazines,  a  nominal 
value  (fixed)  was  chosen  tor  carbon  dioxide,  and  the  value  tor  water  was  computed 
from  the  atmospheric  conditions  as  given  in  Appendix  A. 

Overall,  the  coaputat ion.t  1  procedure  requires  a  pool  temperature  estimate. 

In  all  this  work,  the  first  est  im.it  e  was  taken  to  be  the  air  temperature,  Ta.  Based 
upon  this  estimate,  a  set  of  appropriate  physical  properties  was  deduced,  and  an 
evaporative  flux  determined.  Finally,  with  estimated  properties  and  estimated  flux. 
Equation  (5)  could  be  solved  for  the  pool  temperature. 

An  overall  lntcration  fol-lowcd,  based  upon  revised  estimated  pool  temperatures 
until  satisfactory  agreement  existed,  at  which  time  the  corresponding  evaporative 
flux  was  taken  as  the  correct  value,  and  the  total  evaporation  rate  was  determined 
as  the  product  of  the  flux  times  the  pool  area. 

Dispersion  calculations  used  the  evaporation  rate  as  the  source  strength  for 
the  hydrazine  materials.  The  basic  assumptions  were  that  the  source  was  a  ground- 
level  point  source  with  no  plume  rise,  arid  in  case  of  two  toxic  materials  present 
simultaneously  (aerozinc  50),  the  toxicity  was  considered  additive,  but  not  syner¬ 
gistic  after  the  method  of  the  American  Conference  of  Industrial  Hygenlsts  (14). 

Any  reaction  of  the  evaporated  hydrazine  as  derivatives  in  air  was  treated  as 
a  first-order  decay.  The  equations  for  dispersion  are  as  given  by  Turner  (15)  and 
as  used  in  Illc  and  Springer  (}).  The  basic  relation  for  centerline  concentration 
for  any  distance  x  downwind  is: 

c(x,o)  “  q*A-P-  ■  e»P  ”  <krx/u> 

W0„0,u 


(7) 
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where: 

k  -  atmospheric  decay  constant 
Ap  •  pool  area 

Oy  *  crosswind  diffusion  coefficient 
0Z  •  vertical  diffusion  coefficient 


At  any  distance  y  away  from  the  centerline,  the  concentration  is  given  by: 


c(*.y>  "  cu.o)  e*1’  -  l/2(iy)2 


(7, a) 


The  diffusion  coefficients,  Oy  and  af  are  functions  of  the  downwind  distance, 
x.  The  values  used  In  tills  work  are  those  given  In  llle  and  Springer  (3)  and  are 
as  recommended  by  the  Env  1  r onsient  a  1  Tt  t  e  t  i  or.  \  ,*m  v  . 

For  the  case  >f  more  than  on*  lost.  mater  a.,  th«  existence  of  a  hazardous 
condition  may  be  determine.:  i.  follow  >r  tmi.  .  ‘ti'mentc: 


%  ;  cj 

‘  i-i  Li 

Ct 

where  the  ratio  "  osbin.  a.-.u 

Li  ”  Exposure  limit 


.exponent  1 


C|  ”  Cone ent t at lo i  of  component  1 

If  the  value  of  ts  less  than  unity,  tin-  exposure  limit  is  not  exceeded. 
“T 

In  use,  the  dispersion  calculation  was  carried  out  by  evaluating  F.quatlons 


(7)  and  (7, a)  at  particular  sets  of  ■  and  y  values  for  each  toxic  component,  then 
evaluating  liquation  (8)  to  determine  If  a  hazard  existed.  Practically,  the  down¬ 
wind  distance  was  successively  incremented,  and  at  each,  the  crosswind  distance  was 
incremented  until  the  combined  hazard  ratio  was  less  than  unity.  The  computation 
was  completed  when  the  centerline  (y  *  0)  concentrat ions  (with  successive  downwind 
distances)  showed  a  combined  hazard  ratio  of  less  than  unity. 
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Computational  Results 

The  desired  output  of  the  computing  program  was  an  identification  of  the  hazard 
corridor  for  a  particular  spill.  It  was  considered  that  the  hazard  corridor  was 
defined  if  the  maximum  downwind  distance  and  the  maximum  crosswind  distance  of  the 
hazard  area  were  determined. 

The  computed  output  could  have  been  used  to  determine  the  shape  of  the  hazard 
area,  but  there  are  such  a  large  number  of  variables  In  the  conditions  that  it 
seemed  appropriate  to  simply  define  a  length  (maximum  downwind  distance  of  hazard) 
and  width  (maximum  crosswind  distance  of  hazard). 

The  several  variables,  atmospheric  conditions  as  well  as  spill  dimensions  were 
arbitrarily  set  for  a  "base  case"  for  each  of  the  propel la~t  mixtures  considered. 
Then,  in  order  to  determine  the  effects  of  the  several  variables  independently,  a 
aeries  of  computations  were  made  wherein  the  selected  variable  was  varied  incremen¬ 
tally  from  the  base  case,  and  the  results  determined  in  the  form  of  "correction 
factors".  The  correction  factors  derived  from  these  calculations  arc  plotted  in 
the  figures  following.  Tabulated  values  arc  presented  in  Appendix  F.. 

Table  I  defines  the  "base  case"  conditions,  and  Table  II  shows  the  results  of 
the  base  case  calculation. 
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DEFINITION  OF  hAbl  CAlE  SPILL  CONDITIONS 
Spill  size:  16  ■'  (csrlMd) 

Spill  area :  ISO  m* 

Spill  diameter*:  20  m 

Air  temperature:  25*C 

Ground  temperature:  25*0 
Wind  speed:  3  m/a  (6.7  miles/hr) 

Stability  class:  C 

Relative  humidity:  60* 

Latitude:  30*  North**,  Longitude:  100*  West** 

Date:  266**  (September  23),  Time:  10:00  a.m.** 

Cloud  cover:  30X** 

•Spill  "diameter"  is  actually  the  distance  across  the  pool 
in  the  direction  of  the  wind. 

**Thc  date,  time,  location  and  cloud  cover  conditions  result 
in  an  insolation  (solar  flux)  of  440  w/m^. 
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BASE  CASE  DISPERSION  RESULTS 


Mixture 

Mass 

of 

Spill 

(**) 

Max lmum 

Downwind 
Distance 
of  Hazard 

Maximum 

Crosswind 
Distance 
of  Hazard 
(m) 

Amount  of 
Water  to 
Eliminate 
Hazard* 
(Rg) 

Hydrazine  :  water 

36070 

360 

35 

216000 

(70  :  30) 

Hydrazine  :  I’DMH 

32620 

630 

58 

780000 

(50  :  50) 

HNH  "pure" 

31330 

1630 

138 

4480000 

•The  uount  of  water  which  would  be  required  to  sufficiently  dilute  the  spill 
ao  that  It  would  no  longer  present  a  hazard  was  computed  by  assuming  a  hazard 
would  no  longer  exist  If  the  downwind  distance  of  the  hazard  were  no  more  than 
the  pool  diameter. 


Since  neither  the  downwind  nor  crosswind  distances  arc  linear  in  pool  area, 
plots  of  those  distances  as  functions  of  pool  area  arc  shown  in  Figures  1  and  2. 
The  figures  3  through  10  arc  correction  factors  to  the  data  presented  In  Figures 
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In  order  to  provide  a  method  of  rapid  field  estimation  ol  the  extent  of  a 
hazard  corridor,  the  maximum  distances  of  hazard  presented  in  Figures  1  and  2 
may  be  multiplied  by  one  or  more  of  the  correction  factors  presented  in  the 
following  Figures  3  through  10. 

Comparisons  between  the  computed  result  and  the  estimated  result  (using 
the  correction  factors)  have  indicated  that  the  maximum  difference  between  com¬ 
puted  and  estimated  results  would  probably  not  exceed  20  percent,  although  all 
possibilities  have  not  been  explored.  Also,  for  small  spills,  the  use  of  dis¬ 
tance  increments  of  ten  meters  may  Introduce  some  additional  rounding  and  trun¬ 
cating  errors.  An  error  of  25  percent,  if  that  were  the  case,  would  not  be  an 
excessive  error  for  such  a  computation,  and  could  actually  be  safely  ignored 
when  it  is  compared  to  the  possible  error  in  the  dispersion  calculations. 

Turner  (15)  points  out  that  the  dispersion  calculation  result  will  be  "within 
a  factor  of  three"  of  the  actual  concentration  when  the  terrain  and  circumstances 
agree  reasonably  well  with  the  assumptions  of  the  model. 

An  example  of  the  use  of  the  precomputed  results  follows: 

Example  case : 

Spill  of  250  kg  of  MMH  with  a  spill  area  of  10  m* 

Rci.  Humidity  «  402 
Cloud  Coverage  -  02 
Wind  Speed  -  2  m/s 
Stability  Class  ■  B 
Air  Temperature  ”  30*C 
All  other  conditions  similar  to  "Base  Case". 

From  Fig.  1,  the  downwind  maximum  Is  250  m 


From  Fig  2,  the  crosswlnd  maximum  Is  21  m 


Air  temperature  correction.  Fit;.  3,  1*  1.03 

Stability  class  correction:  Flu.  6,  Downwind,  O.b;  Fig.  7,  Crosswind,  0.97 
Wind  velocity  correction.  Kit;.  <*.  is  I*®4* 

Cloud  cover  correction.  Fig.  9,  Is  1.00 
Relative  humidity  correction.  Fig.  10,  Is  0.98 

Estimated  values: 

Downwind  max.  -  250  x  1.03  x  0.6  x  1.08  x  1.00  x  0.98  •  163  m 
Crosswind  max."  27  x  1.03  x  0.97  x  1.08  x  1.00  x  0.98  -  29  ns 

The  computed  values  for  this  case  are: 

Downwind  max.  *  170  cs 
Crosswind  max.  ■  28  a 
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Limitations  of  the  Model: 

While  the  model  elements  which  have  been  assembled  represent  the  state- of- 
the  art  In  predicting  evaporation,  there  are  several  areas  ol  weakness  among 
the  components  of  the  model.  Each  of  these  weaknesses  are  considered  below, 
although  in  general  it  is  not  possible  to  evaluate  the  magnitude  of  error  which 
would  follow  from  possible  Inapplicability  of  assumptions  or  models. 

Dispersion  Estimates 

As  was  mentioned  In  the  proceeding,  the  dispersion  estimation  may  have  errors 
of  substantial  magnitude.  Turner  (15)  suggests  for  level  terrain,  and  for  center¬ 
line  concentrations  close  to  the  source,  that  the  estimated  value  is  within  a 
factor  of  three,  which  seems  to  indicate  a  possible  error  of  some  200  percent. 
Compared  to  this,  the  other  errors  in  the  model  are  not  very  significant. 

Hass  Transfer  Model 

The  mass  transfer  model.  Equation  (1),  taken  from  HacKay  (10)  is  intended 
as  a  coefficient  for  a  pure  material  evaporating,  or  if  not  pure,  a  well  stirred 
mixture  in  the  liquid  pool.  The  reason  is  that  the  model  assumes  all  the  resis¬ 
tance  of  mass  transfer  is  in  the  vapor  film.  It  is  probable,  however,  that  a 
considerable  liquid  phase  resistance  will  be  present  in  most  of  the  practical 
applications  considered  in  this  study;  even,  perhaps,  the  pure  MMH,  since  it  is 
hygroscopic  and  would  develop  a  water  rich  layer  immediately  adjacent  to  the 
interface  as  water  is  absorbed  from  the  atmosphere. 

Although  it  would  be  difficult  to  quantify  the  effect  of  the  error  in  this 
instance,  it  would  be  most  significant  for  components  which  evaporate  rapidly. 

Equilibrium  Concentrations 

The  binary  and  ternary  system  vapor-liquid  equilibrium  data  are  necessary 
for  predicting  the  concentration  driving  force.  All  the  vapor-llquld  equilibrium 
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data  which  were  available  lor  this  study  were  at  one  atmosphere  total  pressure 
ol  transported  component,  that  Is,  at  the  bubble  point  (boiling  point)  of  the 
materials.  Actually,  the  partial  pressures,  and  hence  the  temperatures,  will 
be  much  lower  than  that  of  the  data. 

To  use  this  data.  It  was  assumed  that  the  activity  coefficients  were  func¬ 
tions  only  of  composition.  For  the  binary  system  of  hydrazine  and  water,  the 
relation  of  activity  coefficient  to  composition  could  be  treated  as  a  polyno¬ 
mial.  However,  for  the  ternary  system  of  hydras lne-UDMll-wa ter ,  the  binary 
data  had  to  bo  fitted  to  an  integrated  form  of  the  Cibbs-Duhem  Equation,  in  this 
case,  the  Harqules  Equation  (7,  8)  In  order  to  use  the  ternary  form.  The  binary 
data  did  not  fit  the  Harqules  Equation  well,  and  so  there  must  be  some  error 
In  the  ternary  predictions. 

In  the  case  of  the  MMH-water  system,  there  were  no  data  available  at  all, 
and  were  estimated  as  being  intermediate  between  the  other  two  binaries,  hydra¬ 
zine-water  and  UDMli-watcr. 

The  error  from  these  assumptions,  especially  for  the  system  MMH-water  could 
be  as  high  as  seventy-five  percent. 

Other  Errors 

Most  of  the  remaining  errors  are  in  the  properties  evaluations,  such  as 
viscosity,  density,  thermal  conductivity  and  so  forth  of  the  mixtures.  These 
errors  will  be  relatively  insignificant  because  of  the  much  larger  uncertainties 
of  the  dispersion  model,  and  perhaps,  the  equilibrium  concentrations. 

Generally  these  properties  were  computed  by  recognized  correlations,  and 
will  give  results  with  less  than  ten  percent  error. 
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APPENDIX  A 

Property  Evaluations  and  Miscellaneous  Computational  Methods 
I.  Properties  of  Liquid  Mixtures: 

The  liquid  properties  of  interest  in  these  computations  were: 

Density,  p 
Viscosity,  u 
Thermal  Conductivity,  k 
Heat  Capacity,  Cp 

Coefficient  of  Thermal  Expansion,  8 

The  methods  of  computation,  including  the  effect  of  temperature  on  the  property 
were  generally  as  recommended  by  Reid  and  Sherwood  (16). 

Density 

For  a  liquid  mixture*  the  volume  was  assumed  to  be  additive,  hence  the 

density  was  the  mass  fraction  weighted  sum,  thus: 
n 

Pm  ■  £xiPi  (A-l) 

1-1 

where:  X*  “  mass  fraction  of  component  1  and 

Pi  -  density  of  component  1. 

Individual  expressions  for  the  density  of  each  component  as  a  function  of 
temperature  were  obtained  as  polynomials,  either  as  presented  directly  or  by 
reduction  of  published  data. 

For  hydrazine  (9) 

p  -  1210.78  -  0.62668 (T)  -  4. 5284 x  10‘4(T2)  (A-l-a) 

For  UDMH  (9) 

P  -  1060.41  -  0.77507 (T)  -  4.8648 x  10_4(T2)  (A-l-b) 

For  ttOi  (9) 


p  -  1150.34  -  0.93949(T) 


(A-l-c ) 
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For  water  (17) 

P  -  1054. A  -  0. 19241 (T )  (A-l-d) 

For  the  above  expressions ,  the  density  p  is  in  kg/m-*  and  the  temperature 
ia  in  kelvins. 

Viscosity 

The  liquid  mixture  viscosities  were  determined  from  methods  presented  by 
Reid  and  Sherwood  (b). 

log[log(10yB)  ]  -  -  2.900  (A—  2 ) 

In  this  expression,  um  -  mixture  viscosity  in  centlpoise 
Pb  ■  mixture  density  in  g/cm* 
x*  *  mole  fraction  of  component  1 
Hi  -  molecular  weight  of  component  i 
I|  «  viscosity  constant  for  component  1 
for  hydra  tine,  l  -  84.8 

for  MMH,  1  -  140.4 

for  UDMH,  I  -  196.0 

for  water,  1  -  59.8 

The  temperature  effect  on  viscosity  was  treated  by  the  "Souder ' a  Method" 
as  given  in  Reid  and  Sherwood  (16)  and  this  is  Included  in  the  foregoing  equa¬ 
tion  for  the  viscosity  of  liquid  mixtures. 

Thermal  Conductivity 

Thermal  Conductivity  of  the  liquid  mixtures  was  treated  bv  mole  fraction 
weighting,  thus: 

m 

k.  -  IkiXi  (A-J) 

1*1 

where:  k^  -  thermal  conductivity  of  component  1 
X|  *  mole  fraction  of  component  1 
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The  individual  component  conductivities  were  determined  as  temperature 
function  polynomials,  either  as  directly  presented  or  by  reduction  of  published 
data  (9.  17). 

For  hydrazine :  k  -  0.088142  +  2. 742 x 10"3(T)  -  4.677  x  10_6(T2)  (A-3-a) 

For  UDMH:  k  -  0.28841  -  4.262 x  10~4(T)  (A-3-b) 

For  MMH:  k  -  0.14236  +  9. 2048 x  10‘4(T)  -  1 . 9014 x  10~6(T2)  (A-3-c) 

For  water:  k  -  0.3427  +  9.178 x  10_4(T)  (A-3-d) 

The  thermal  conductivities  are  in  w/mK  when  the  temperature,  T,  is  in 
Kelvins. 


Heat  Capacity 

The  heat  capacity  values  of  the  liquid  mixtures  were  computed  as  mass 

weighted  averages,  thus: 

m 

CD  "  £CPi*i  <A-*) 

“  1-1 

where:  Cp  -  heat  capacity  of  component  1 
Xj  -  mass  fraction  of  component  i 

The  component  heat  capacities  were  computed  from  temperature  polynomials, 
either  as  directly  presented  or  as  deduced  by  data  reduction  (9,  17) 


For 

hydrazine: 

CP 

1234.8 

+  8. 4488 (T)  -  7.7651  x  10~3(T2) 

(A-4-a) 

For 

UDMH: 

CP  " 

1703.3 

+  3.6978(T) 

(A-4-b) 

For 

MMH: 

CP“ 

2731.3 

-  0. 07 2316 (T)  ♦  2.4745  x  10‘3(T2) 

(A-4-c ) 

For 

Water: 

CP  " 

2512  ♦ 

6. 147 (T) 

(A-4-d) 

where:  Cp  -  heat  capacity,  J/kg*K 
T  -  temperature,  K 

Coefficient  of  Thermal  Expansion 

The  coefficient  of  thermal  expansion  was  computed  by  noting  the  change 
in  density  with  temperature: 
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B - £L.-.BL _ 

2(TX  -  T2)PiP2 


(A—  5 ) 


where:  Pj,  P2  “  liquid  mixture  density  at  temperatures  of  Tj  and  T 2  respectively. 

Thla  function  was  computed  each  time  needed  by  finding  the  density  at  two 


different  temperatures  in  the  temperature  region  of  interest. 


II.  Properties  of  Vapor  Mixtures: 

The  vapor  properties  necessary  for  computation  of  the  properties  of  the 
alr/llquid  surface  film  included: 
d  ensity 
viscosity 
heat  capacity 
thermal  conductivity 
dlf  f uslvlty 
Schmidt  no. 

For  the  hydrazines,  the  basic  Information  on  vapor  phase  viscosity,  heat 
capacity  and  thermal  conductivity  was  taken  as  polynomials  from  the  source  pre¬ 
viously  cited  (9).  For  water,  these  properties  were  generally  taken  from  Kern 
(17).  Additionally,  the  vapor  film  properties  are  predoml natelv  those  of  air. 
Air  properties  were  taken  generally  from  Perry  (12). 

The  density  was  computed  by  assuming  the  validity  of  the  Ideal  gas  law. 


Viscosity 


J1 


For  water: 
For  air: 


Hv 

Uv 


33.333  x  10'8(T)  -  1 x  10“7 

1.45  x  10~fe(T)1-5 
(T  ♦  116) 


where:  Uv.l  “  vapor  phase  viscosity  in  kg/m- 5 

y^  •  mole  fraction  of  component  i  in  vapor 
“  molecular  weight  of  component  1 
T  •  temperature  in  Kelvins 


(A-6-d ) 
(A-6-e) 

(A-7) 


Heat  Capacity 

Vapor  phase  heat  capacities  were  computed  by  polynomial  expressions  of 
the  form  Cp^  “  a^  ♦  bjT  +  CjT-.  For  the  components,  the  following  constants 

were  used: 


Component 

a 

b 

£ 

hydrazine 

1493.7 

3.3133 

1.0209  x  10" 3 

UDMH 

264.43 

4.548 

-1.5146  x  10~3 

MMH 

355.22 

4.1338 

-1.34725  x  10" 

water 

1656.9 

0.61297 

1.0657  x  10“  3 

air 

970.69 

0.067864 

1.6569  x  10‘3 

The  heat  capacity  of  the  mixture, then  was: 


m 

Cp^  *  L  yiCpi  (A— 8) 


The  above  constants  give  heat  capacity  values  in  J/kg-K  when  temperatures  are 


in  Kelvins 
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Thermal  Conductivity 

The  vapor  phase  thermal  conductivity  was  estimated  by  methods  recommended 
by  Reid  and  Sherwood  (16),  as  viscosity  and  heat  capacity  functions.  The  com¬ 
ponent  thermal  conductivities  were  computed  by  the  Eucken  approximation  as  given 
by  Perry  (12). 

ki  "  MCpi  +  10393/Mi)  (A-9) 

where  is  the  component  molecular  weight.  For  the  mixture,  the  thermal  con¬ 
ductivity  was  computed  by: 


5  .  m1/3 

t  yikiMi 

v  «l/3 

l  y^ 

i-1 


(A-10) 


Dlffusivlty  and  Schmidt  Nvimber 

The  values  for  dlffusivlty  of  each  component  in  air,  and  the  corresponding 
Schmidt  number  were  computed  by  mt thods  given  by  Sherwood,  Pigford  and  Wilke 
(18).  These  were  precalculated  as  temperature  functions  for  all  the  components 
of  interest,  so  In  use,  the  component  and  the  temperature  was  all  that  was  neces 
sary  to  compute  a  value  for  Schmidt  number  or  the  binary  dlffusivlty.  No  allow¬ 
ance  was  made  for  the  mixture  having  more  than  two  components.  Each  dlffusivlty 
each  Schmidt  number,  and  hence,  each  mass  transfer  coefficient  was  determined 
as  if  it  were  the  only  component  being  transfered. 


III.  Computation  of  Mass  Transfer  Flux  in  the  Absence  of  Wind. 

In  order  to  account  for  a  small  mass  transfer  effect  that  would  occur  in 
the  absence  of  wind.  Equation  (1)  had  to  be  modified  since  it  can,  as  written, 
account  only  for  the  turbulent  transport. 

Sherwood,  et.al.  (18)  present  a  relation  for  predicting  flux  at  the  sur¬ 
face  of  a  sphere  under  molecular  diffusion.  This  relation  can  be  rearranged  to 
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give  a  static  mass  transfer  coefficient,  k8  of: 
DP 


U  m 

“  RTrP 


BM 


(A- 11) 


where:  D  -  diffusivity, 

P  "  pressure, 

R  -  gas  constant, 

T  -  temperature ,  and 
r  -  radius  of  the  sphere. 

Pgjj  “  mean  effective  pressure  of  inert 
The  pool  surface  was  likened  to  the  surface  of  a  sphere,  and  the  equivalent 
radius,  in  terms  of  the  pool  area  was  found.  The  resulting  expression  is: 

k“  ■p'hh*^  lA_12> 

where  Ap  is  pool  surface  area.  This  expression  for  static  transfer  was  added 
to  the  expression  for  turbulent  transport  given  in  Equation  (1). 

IV.  Computation  of  Equilibrium  Vapor  Compositions 

The  basic  equation  used  to  compute  vapor  equilibrium  composition  was: 

Yi  -  (A-l  3) 

where:  y^  ■  vapor  phase  mole  fraction  of  component  1 

«  activity  cocfficitnt  of  component  1 
Xj  ■  Required  phase  mole  fraction  of  component  i 
PV>1  ”  vapor  pressure  of  component  1 
P  »  total  pressure 

The  activity  coefficient  is  a  function  of  the  liquid  phase  composition. 

The  equation  used  in  this  work  was  the  three-suffix  Marqules  equation.  For  a 
binary  mixture: 


log  ■  X21  Aj2  ~  2X^21  -  A12)J 
log  I2  "  X^lAji  -  2X2(A|2  “  A2i) ) 


(A-U,a) 

(A-lA.b) 


The  coefficients  for  a  ternary  mixture  can  be  deduced  from  the  binary  data 


if  the  constants  A33»  Aj^,  Ai3>  AJ1>  A23  a,w*  A32  are  known,  in  which  case: 
log  yx  -  I a12  -  2*1<A21  -  a12>J  +  X  j  [  Ai3  +  2Xl(A31  -  Au)l 
*2*3 1 1^2 (A21  ^  a12  *  A31  *  Ai3  "  a23  "  A3->) 

♦  ^i(A21  “  a12  +  a31  ~  a13^  +  "  Xj)(A23  -  A^)] 

For  /g*  £he  subscripts  are  rotated,  1  becoming  2,  2  becoming  3  and  3  becoming  1. 
For  jTj,  the  subscripts  are  rotated  once  more. 


p 

BM 

Pv 

Pr 

*a 

^n 

«lr 

Qm 

r 

R 

Sc 

T 

Ta 

Tg 

Tp 

u 

X 

x 

y 


pressure 

mean  pressure  of  non-diffusing  components 
vapor  pressure 
Prandtl  No. 

heat  flux  from  atmospheric  radiation 

heat  flux  from  convective  exchange  between  pool  and  gro<  :id 

heat  flux  from  convective  exchange  between  pool  and  air 

latent  heat  of  vaporizing  or  condensing  materials,  expressed  as 

reradiation  heat  flux  from  pool 

solar  flux 

evaporative  t lux 

radius  of  sphere  equir  to  pool 

universal  gas  constant 

Schmidt  number 

temperature 

temperature  of  the  air 
temperature  of  the  ground 
temperature  of  the  pool 
wend  velocity 

distance  across  pool,  also  mole  fraction  in  liquid 
mass  fraction  In  liquid 

mole  fracti«-  vapor,  also  distance  away  from  centerline 


6 

* 

V 

P 

a 


coefficient  of  expansion 
activity  coefficient 
viscosity 
density 

Stef  f an-Haltzman 
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APPENDIX  C 

Description  ot  the  Computer  Program 

A  computer  program  was  written  (see  Appendix  D)  to  estimate  the  evaporation  flux 
and  the  ground  level  concentrat ions  which  result  from  the  dispersion.  Figure  11 
shows  a  simplified  flow  diagram. 

The  essential  feature  of  the  model  is  that  when  the  temperature  of  the  pool 
has  been  satisfactorily  estimated,  the  flux  can  be  readily  deduced.  However,  since 
the  evaporation  flux  has  a  significant  bearing  on  the  pool  temperature,  an  itera¬ 
tive  solution  is  indicated. 

The  program  makes  an  initial  assumption  that  the  pool  temperature  is  the  same 
as  the  air  temperature.  After  this  assumption,  all  the  properties  and  coefficients 
are  evaluated  and  the  evaporative  flux  computed.  When  the  flux  is  computed,  the 
energy  balance  is  solved  for  the  pool  temperature  usinR  the  Newt on-Raphson  method. 

If  the  new  solution  for  the  pool  temperature  is  used  isEsediately  as  a  revised  esti¬ 
mate,  then  under  some  conditions  the  overall  s olution  diverge*.  Therefore,  a 
damping  factor  is  used  to  obtain  a  new  temperature  assumption. 

When  the  solution  has  converged  satisfactorily,  then  the  last  value  of  the 
evaporative  flux  is  taken  as  the  correct  value  and  the  program  proceeds  to  com¬ 
pute  an  evaporation  rate. 

The  evaporation  rate  is  assumed  constant  for  one  second,  and  the  amount  of 
each  component  transferlng  in  the  one  second  period  is  subtracted  from  the  pool 
inventory,  so  a  new  pool  composition  may  be  computed. 

After  each  one  second  of  evaporation,  a  new  composition  is  used.  The  compu¬ 
tation  proceeds  by  su«lng  the  one  second  amounts  for  16  seconds  (0.01  hour).  The 
average  rate  over  the  36  seconds  Is  then  used  for  the  dispersion  calculations  as 
the  source  strength. 

For  each  ten  meters  down  wind,  the  dispersion  subroutine  computes  the  ground- 
level  concentration  at  the  centerline  and  at  Incremental  distances  away  from  the 
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centerline  until  the  result  shows  the  concentration  to  be  below  the  hazard  level. 
That  is,  the  downwind  distance  is  Incremented  by  ten  meters,  then  the  crosswind 
distance  is  incremented  and  so  forth,  until  finally  the  centerline  concentration 
is  found  to  be  below  the  hazard  level. 

Various  output  information  is  available;  the  actual  output  can  be  whatever  is 
desired . 

Among  the  variables  computed  which  would  be  of  interest  as  output  are: 

CrosswinJ  distance  of  the  hazard  as  a  function  ol  downwind  distance; 

Maximum  crosswind  distance  of  the  hazard; 

Maximum  downwind  distance  of  the  hazard; 

Evaporation  rate  of  each  component; 

Liquid  pool  composition. 

All  of  the  above  are  computed  for  each  0.01  hour 

The  input  data  required  are: 

The  composition  of  the  original  spill,  mass  fraction  of  hydrazine,  UDMH, 
water  and  MM11; 

The  mass  (kg)  of  the  spill,  the  "diameter"  (downwind  distance  across 

2 

the  pool)  of  the  pool  (m)  and  the  pool  area  (m‘); 


The  time  zone  (from  Greenwich),  latitude,  longitude  (degrees),  date, 
(Julian  day),  time  (24  hr,  decimal  hours),  the  air  temperature  (*C), 
ground  temperature  (*C),  wind  velocity  at  10  m  elevation  (m/s),  rela¬ 
tive  humidity  (percent),  cloud  coverage  (percent),  and  the  stability 
class  integer,  (1,  2,  3,  4,  3,  b  for  A,  B,  C,  0,  E,  F,  respectively). 

Since  the  program  computes  the  evaporation  and  dispersion  for  one  hour,  new 
weather  data  may  be  read  in  for  succeeding  hours  up  to  24. 
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Sub rout lnos 

The  main  program  calls  the  following  subroutines  as  required: 

Subroutine  S01LM:  Returns  ground-to-pool  heat  transfer  coefficient  as  a  function 
of  the  pool  composition,  amount,  area,  temperature  and  the  soil  temper¬ 
ature. 

Subroutine  SKY:  Returns  long  wave  length  radiation  flux  and  the  partial  pressure 

of  water  vapor  from  the  atmosphere  as  functions  of  the  air  temperature, 
relative  humidity,  and  cloud  cover. 

Subroutine  SOLAR:  Returns  solar  flux  intensity  as  a  function  of  date,  time,  lati¬ 
tude,  longitude  and  cloud  cover. 

Subroutine  AJUUi:  Returns  air  to  pool  heat  transfer  coefficient  as  a  function  of 
mass  transfer  coefficients,  Sctimidt  numbers,  vapor  composition,  air 
temperature  and  pool  temperature. 

Subroutine  SCMIT :  Returns  Schmidt  numbers  and  d if f uslvlt les  for  the  transfcrcd 
components  in  air  as  functions  of  the  temperature  of  the  pool  and  of 
the  air. 

Subroutine  GAMMA :  Returns  equilibrium  partial  pressures  for  all  transfered  components 
as  functions  of  composition  and  temperature  of  the  pool. 

Subroutine  D1SPER:  Computes  and  prints  the  results  of  dispersion  calculations. 

The  Input  is  the  evaporation  rate,  which  is  the  source  strength,  and 
weather  data. 

Appendix  D  contains  a  listing  of  the  program  and  a  sample  of  the  output. 


Variable  Name 


Variable  Description 


Xll,  XI 2,  XI],  X14 

HASS 

DLA 

AREA 

DATE 

TIME 

TA 

TG 

U 

RHUM 

XLAS 

mprijtt 

DELT 

DUS 

CN 

POO  LAB 

PT 

TIN 

XI.  X2.  X).  X4 

TP  V 

TCC 


Initial  composition  mass  2  of  hydrazine,  UDMH,  water, 
and  MMH,  respectively 

mass  of  the  propellant  spill,  kg 

diameter  of  the  propellant  spill,  meters  (downwind 
distance) 

area  of  the  spill,  m" 

Julian  day  (1  to  36b) 

time,  hours  after  midnight  (Oto  24) 

air  temperature,  *C  (sometimes  K,  Internally) 

ground  temperature,  *C  (sometimes  K,  Internally) 

wind  speed,  m/s 

relative  V-aid lty ,  percent 

atmospheric  stability  class,  numerical,  1  "  A,  2  •  B, 
etc . 

output  printing  scheme  variable 

time  constant  In  seconds  used  in  calculations 

pool  emlsslvlty 

ground  roughness  factor  and  temperature  factor  In  the 
atmosphere 

pool  absorptivity 

atmospheric  pressure.  Pa 

initial  time  in  hours  after  midnight 

mass  fraction  of  hydrazine,  UDMH,  water,  and  MMH, 
respectively 

computed  temperature  of  the  liquid  pool,  K 

turbulent  contribution  to  the  mass  transfer  coefficient. 


Variables  Used  In  Computer  Program  (continued) 


Variable  None 

PBM1,  PBM2,  PBM3, 
PBM4 ,  PBM5 

TCI,  TC2,  TC3, 

TC4,  TC5 

HVAP1,  UVAP2 ,  HVAP3, 
HVAP4 

TPT 

EVAP1,  EVAP2,  EVAF3 , 
EVAP4 

C02 

qvAP 

CONST 

XM1 ,  XM2,  XM4 

XJ1,  XJ 2 ,  XJ3,  XJ4 

SUM 1 ,  SUM2,  SUM3, 
SUM4 

AVC1,  AVC2 ,  AVC3, 
AVG4 

XM 

YM 

CCE 

CC1,  CC2,  CC4 
CCY1 ,  CCY2,  CCY4 


Variable  Description 


mean  partial  pressures  of  non-diffusing  material.  Pa 


mass  transfer  coefficient  for  hydrazine,  UDMH,  water 
MMH,  and  CO2,  respectively 

heat  of  vaporization  for  hydrazine,  UDMH,  water,  and 
MMH,  respectively,  J/kg. 

temporary  pool  temperature  calculated  by  using  Kaphson 
Newton  method 

evaporation  flux  of  hydrazine,  UDMN,  water,  and  MMH, 
respectively 

Diffusing  flux  of  CO2  into  the  liquid  pool,  kg  mole/m* 

evaporative  heat  loss  of  the  propellant  pool,  J/m-  s 

solar  flux  +  radiative  heat  gain  from  atmosphere  - 
evaporative  heat  loss,  W/nr 

■ole  fraction  In  propellant  spill  of  hydrazine,  UDMH, 
and  MMH,  respectively 

weight  percent  of  hydrazine,  UDMH,  water,  and  MMH, 
respect ively 

algebraic  sum  of  the  evaporation  rates  of  the  respec¬ 
tive  components  over  a  36  second  time  interval,  kg/s 

average  evaporation  rate  over  a  36  second  time  inter¬ 
val  for  hydrazine,  UDMH,  water,  and  MMH,  respectively, 

kg/s 

downwind  distance,  m 
crosswind  distance,  m 

centerline  concentration  equivalent  (dimensionless 
ratio) 

centerline  concentration  in  wg/m^ 
concentration  at  YM,  wg/*} 


XMDEL 


Increment  of  downwind  distance,  m 


Variable  Name 


Variable  Descr ipt ion 


RXK1 ,  RXK2,  RXXh 

Ql.  Q2.  Q4 

XX 

S1CZ 

SICY 

HAST 

I  NCR 

IF 

SCI,  SC2,  SC  3  ,  SC  4 ,  $C5 

D1FF1,  D1KF2,  DIFF3 , 
DIFF4,  D IFF 5 

GAM ( 1 ) ,  CAM(2), 

GAM( 3) ,  GAM ( 4 ) 

VP1,  VP 2,  VP3 ,  VP4 
Ml,  W2,  M3,  M4 

ac 

Q 

T 

DN1 ,  DN2,  DN3,  DN4 

DNLM 

DNMTG ,  DNMTP 

BETA 


reaction  rate  constants  (first  order),  sec 
emission  rates  of  components,  kg/s 
downwind  distance,  km 

dispersion  coefficient  in  vertical  direction,  m 
dispersion  coefficient  in  horizontal  direction,  m 
Counting  Integer,  dimensionless 
print  frequency  increment 
average  film  temperature,  K 

Schmidt  number  for  hydrazine,  UDMH,  water,  MMH,  and 
carbon  dioxide,  respectively 

Dlffusivity  of  hydrazine,  UDMH,  water,  MMH,  and  CC>2 
respectively,  in  air,  m“-/s 

activity  coefficient  for  hydrazine,  UDMH,  water,  and 
MMH ,  respectively 

pure  component  vapor  pressure  for  hydrazine,  UDMH, 
water,  and  MMH,  respectively.  Pa 

mass  of  hydrazine,  UDMH,  water,  and  MMH,  respectively, 
*4 

heat  transfer  coefficient  between  the  ground  and  the 
liquid  pool,  W/m*-  K 

total  mass  of  the  spill,  kg 

temperature  used  in  temp  dependent  expressions,  K 

liquid  density  of  hydrazine,  UDMH,  water,  and  MMH, 
respectively,  kg/m^ 

density  of  the  liquid  mixture,  kg/m^ 

liquid  density  of  the  mixture  at  the  ground  and  pool 
temperature,  respectively,  kg/mJ 

expansion  coefficient,  K  * 


Variable  Name 


Variable  Description 


film  temperature  at  the  soil-pool  interface,  K 
density  of  the  liquid  mixture  at  the  film  temp,  kg/m^ 


TK1,  TK.2,  TKJ ,  TK4 


thermal  conductivity  of  liquid  hydrazine,  UDMH,  water 
and  MMH,  respectively,  W/m  K 


CPI,  CP2,  CP3 ,  CP4 


total  number  of  moles  in  the  liquid  pool,  kg  mole 

liquid  heat  capacity  of  hydrazine,  UDMH,  water,  and 
MMH,  respectively,  in  J/kg  K 

thermal  conductivity  of  the  liquid  mixture,  W/m  K 
heat  capacity  for  the  liquid  mixture,  J/kg  K 
viscosity  of  the  liquid  mixture,  kg/  m  s 
Prandtl  number  for  the  liquid  mixture 
pool  depth  in  meters 

temperature  difference  between  the  pool  and  ground,  K 
Crashoff  number 
Rayleigh  number 


A,  B,  C 


constants  used  in  the  generalized  heat  transfer  co¬ 
efficient  correlation,  also  in  atmospheric  emissivity 
corrclat ion 


CLOUDS 


cloud  coverage  in  percent 


SKYRAD 


partial  pressure  of  water  in  the  atmosphere.  Pa 
radiant  energy  absorbed  by  the  pool,  W/m‘ 


atmospheric  emissivity 

latitude,  degrees  (4  north  and  -  south) 

longitude,  degrees  (+  west  and  -  east) 

time  zone  number  (hours  behind  Greenwich  mean  time) 


intensity  of  the  total  solar  radiation  incident  on 
pool  surface,  W/m^ 


SOLR 


Variables  Used  In  Computer  Program  (continued) 


Variable  Name 
FI,  F2.  F3,  F4 ,  F5 

Z,  W 

SALT 

CCN 

LON 

ID 

SKY 

PI,  P2,  P3,  P4 

U 

WN(1),  WM(2),  WM(3) , 
VM(A>  . 

Y(l),  Y(2) ,  Y(3)  , 

Y(4).  V( 5) 

DNF 

CP(1),  CP(2)  ,  CP(3) , 
CP(4),  CP(5) 

CPM 

VIS(l),  V1S(2) ,  V1S(3) , 
VIS(4),  VIS(5) 

VISM 

TK(1),  TK( 2) ,  TK(3), 
TK(4).  TK(5) 

TKM 

CPF 


Variable  Description 

variables  related  to  solar  radiation  as  functions  of 
the  date 

direction  cosines  of  direct  solar  beam 
solar  altitude,  radians 
cloud  cover  modifier 

2 

intensity  of  direct  normal  solar  radiation,  BTU/hr  ft" 

intensity  of  solar  radiation  Incident  on  pool  surface, 
Btu/hr • f  t ‘ 

•y 

sky  diffuse  radiation  incident  on  surface,  Btu/hr'ft*- 

equilibrium  partial  pressure  of  hydrazine,  UDMH,  water, 
and  MMH,  respectively,  Pa 

heat  transfer  coefficient  between  the  liquid  pool  and 
the  air,  V/a" -K 

molecular  weight  of  hydrazine,  UDMH,  water,  MMH,  and 
air,  respectively 

average  mole  fraction  In  the  pool-air  film  of  hydra¬ 
zine,  UDMH,  water,  MMH,  and  air,  respectively 

molar  density  of  the  vapor  film,  kg'mole/m^ 

gaseous  heat  capacity  of  hydrazine,  UDMH,  water,  MMH, 
and  air,  respectively  In  J/kg*K 

average  molar  heat  capacity  of  the  vapor  film,  J/kg 

mole'K 

gaseous  viscosity  of  hydrazine,  UDMH,  water,  MMH, 
and  air,  respectively,  kg/m*s 

average  viscosity  of  the  vapor  film,  kg/m-s 

gaseous  thermal  conductivity  for  hydrazine,  UDMH, 
water,  MMH,  and  air,  respectively,  W/m’K 

thermal  conductivity  of  the  vapor  film  mixture,  W/m'K 

average  heat  capacity  of  the  vapor  film  mixture,  J/kg‘K 


Variable  Name 


Variable  Description 


TC  average  mass  transfer  coefficient  of  the  vapor  film 

mixture,  m/s 

SC  average  Schmidt  number  of  the  vapor  film  mixture 

PR  Prandtl  number  of  the  vapor  film  mixture  at  the  film 

temp 
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MAXIMUM  CROSSWISD  AND  DOWNWIND  DISTANCES  AS  A  FUNCTION  OF  SPILL  AREA  WITH 


GROUND  TEMPERATURE  -  25*C 


downwind  crosswind 


downwind  crosswind  downwind  crosswind 


.78 

.81 

.89 

.97 

.00 

.14 

.22 

.25 

.28 

.28 


Table  E-2 

SOKRECT ION  FACTORS  FOR  SOLAR  FLUX 
(Base  Case:  SOLR  -  440  w / m - ) 


Solar 


Hydrazine  :  Water 


Hydrazine  :  UDMH 


MMH 

100Z 


953 


Table  E-3 


CORRECTION  FACTORS  FOR  STABILITY  CLASS 


Stability 

Class 

Hydrazine  :  Water 
(70/30) 

Hydrazine  :  t'DMH 
(50/50) 

MMH 

(1002) 

downwind  crosswind 

downwind 

crosswind 

downwind  crosswind 

A 

0.44  1.0 9 

0.43 

1.03 

0.33 

0.91 

B 

0.64  1.06 

0.63 

1.03 

0.58 

0.97 

C 

1.00  1.00 

1.00 

1.00 

1.00 

1.00 

D 

1.58  0.97 

1.79 

1.09 

2.11 

1.20 

E 

1.97  0.89 

2.75 

1. 16 

3.40 

1.33 

F 

3.28  0.91 

4.87 

1.26 

6.13 

1.56 

Table  E-4 


CORRECTION  FACTORS  FOR  WIND  VELOCITY 


Wind 

Speed 

(■/») 

Hydrazine  :  Water 
(70/30) 

Hydrazine  :  UDMH 
(50/50) 

MMli 

(1002) 

downwind 

crosswind 

downwind 

crosswind 

downwind  crosswind 

1 

1.31 

1.29 

1.21 

1.19 

1.28 

1.25 

2 

1.11 

1.11 

1.06 

1.07 

1.07 

1.07 

3 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

4 

0.94 

0.97 

0.95 

0.97 

0.94 

0.96 

5 

0.92 

0.91 

0.90 

0.93 

0.91 

0.92 

6 

0.86 

0.89 

0.87 

0.90 

0.88 

0.90 

77 


Table  t-5 

COKKKtriUN  FACTORS  FOR  AIK  rOtn 


Air 

Temper JCure 


Hydra*  ini'  ; 

_ LLQ/JPJ _ 


Hydra/ln* 

(SO/SO) 


MMH 

(100Z) 


C  C) 

downwind 

crottitwlnd 

ilovnw  lnd 

crosswind 

downwind 

crosbvind 

35 

1.25 

1.20 

1.05 

1.04 

1.13 

1.12 

30 

1.12 

1.09 

1.02 

1.02 

1.03 

1.03 

25 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

20 

0.91 

0.91 

0.98 

0.99 

0.98 

0.98 

15 

0.81 

0.85 

0.97 

0.98 

0.96 

0.96 

10 

0.77 

0.  78 

0.9b 

0.97 

0.94 

0.95 

;h 


Tabic  E-b 

CORRECTION  FACTORS  FOR  GROUND  TEMPERATURE 


Ground 

Temperature 

CC) 

Hydras  Inc  :  Water 
(70/30) 

Hydrazine  :  UDMH 
(50/50) 

MMH 

(1001) 

d ovnw 1 nd 

crosswind 

downwind 

crosswind 

downwind  cros-.wind 

15 

0.97 

1.00 

0.83 

0.84 

0.91  0.93 

20 

1.00 

1.00 

0.90 

0.91 

0.92  0.93 

25 

1.00 

1.00 

1.00 

1.00 

1.00  1.00 

30 

1.03 

1.03 

1.10 

1.10 

1.10  1.10 

35 

1.03 

1.03 

1.21 

1.21 

1.21  1.20 

lublc  E-7 


CORRECT  ION  FACTORS  FOR  CLOUD  COVERAGE 


C  loud 
Coverage 


Hvdruzitu-  :  Water 
(70/ >0) 


Hydrazine  :  UDMIi 
(SO/ 50) _ 


MMM 

(loot) 


U) 

downwind 

crosswind 

downwind 

croaswind 

downwind 

crobswlnd 

20 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

JO 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

40 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

bO 

1.00 

1.00 

1.00 

1.00 

0.99 

1.00 

SO 

0.97 

0.97 

0.9S 

1.00 

0.99 

0.99 

100 

0.92 

0.91 

0.98 

1.00 

0.98 

0.99 

